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SUMMARY 

The blackbody temperature and the albedo of a planet, and 
the variation of both parameters with latitude, longitude, and time, 
are of great value in understanding the climatic and meteorolog- 
ical conditions of the planet. An unchopped radiometer with a wide 
but restricted field of view is capable of such temperature and 
albedo measurements. Coated thermistors mounted in highly 
reflective cones serve as detectors. Their performance as 
sensor elements is analyzed in detail herein to prove the fea- 
sibility of the measurement. The simplicity of the instrumenta- 
tion and the low information bandwidth required make the experi- 
ment equally attractive for earth satellites and space probes. 
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A LOW-RESOLUTION UNCHOPPED 
RADIOMETER FOR SATELLITES 


INTRODUCTION 


Measurement of the optical and thermal properties of the earth and other planets is 
one of the basic experiments that can be performed by means of satellites and space 
probes. Incoming and outgoing radiation determine the energy budget of a planet and are 
responsible for climate and weather in the broad sense. Measurements in narrow spec- 
tral bands can be used to determine atmospheric compositions and temperatures and to 
give information on surface conditions. 

The goal in this field of research could be detailed maps of the planets, showing the 
visual picture, temperatures, and radiation in characteristic bands. All high- resolution 
maps obtainable by photography, television, and spot-scanning techniques require high 
rates of information transmission and, frequently, storage facilities in the probe. For 
many purposes, especially if overall global or planetary studies are involved, only aver- 
age values over rather large areas (e.g., 300 by 300 miles) are of interest. Even in the 
case of earth satellites, where a high- resolution picture is possible, the relatively low 
information rate and the integrating property of the wide-field detector justify simulta- 
neous measurements by devices with low and high resolution. 

The simplest form of such a wide-field radiation measurement is the heat-balance 
experiment suggested by Suomi and Wexler in connection with the International Geophysi- 
cal Year program (Reference 1). In this procedure the temperatures of small spheres, 
coated to discriminate between solar and terrestrial radiation, are sensed by thermistors. 
However, the omnidirectional device has restrictions in accuracy and resolution. It is 
especially desirable in an image-forming experiment to match the field of view of both 
high- and low-resolution instruments. Interesting correlations between cloud cover and 
radiation balance might be expected. 

The main problem in such a simple, unchopped radiometer is the loss in sensitivity 
resulting from the restrictions in the field of view. A technique that overcomes this 
limitation is analyzed in this paper. 
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METHOD OF MEASUREMENT 


In the suggested version a thermistor is mounted in the apex of a highly reflecting 
cone. The detector can be classified as a thermal detector cooled predominantly by 
radiation. Its temperature is governed by the equation of energy balance, 

e t oT t 4 = kjCrT* + k 2 cr T s 4 + k 2 4aT 0 3 (T s - T t ) + k 3 A<* sin /3 . (1) 

(All symbols are defined in Appendix A.) 

Thermal radiation emitted by the thermistor balances radiation from the earth, 
radiation and heat conduction from the satellite, and finally the earth’s reflected solar 
energy. The coefficients k consolidate optical as well as geometrical factors. The two 
unknowns, i.e., the blackbody temperature T e and the albedo of the earth, can be deter- 
mined if a second linearly independent equation is available. 

In the experiment the second equation is provided by a second thermistor with dif- 
ferent optical properties. One detector is coated black; it is equally sensitive to reflected 
sunlight and long-wave terrestrial radiation. The second detector is coated to be highly 
reflective in the visible and near infrared to about 3 microns (Reference 2). The surface 
appears white to the eye even though its emissivity in the infrared is very high. Since 
99.9 percent of the terrestrial radiation is emitted at wavelengths longer than 4 microns, 
both detectors will show the same equilibrium temperature when they face the dark side 
of the earth. On the illuminated side the temperature of the black detector will rise — in 
contrast to that of the white one, which is hardly affected by reflected sunlight. 

Careful measurement of the temperatures of both detectors and accurate calibration 
of all k values permits the determination of albedo and terrestrial radiation.* It has 
been tacitly assumed that the influence of the satellite temperature T s can be kept small 
(k 2 and k 2 « kj), or the thermistor would read essentially the satellite temperature. 
High thermal isolation of the detectors from the satellite structure is essential to the 
experiment, and the major part of this paper is devoted to that problem. 


GEOMETRY OF THE CONE 


The field of view of a thermistor mounted in a cone can best be visualized by utiliz- 
ing the mirror image of the detector (Figure 1). The conical detector is equivalent to a 


*This is exactly true only if radiation from the earth follows Lambert's cosine law; but 
since the atmosphere, with its peculiar absorption bands in the infrared and its tem- 
perature gradients, does not emit as a diffuse surface, certain corrections have to be 
made in the interpretation of data. 
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Figure 1 - Equivalent detector, imaging of thermistor 

spherical one that is restricted in its field of view by an aperture equal to the base of the 
cone. Reflection losses could be taken into account by assuming a reduced emissivity for 
each image. It is justifiable to consider a sphere of uniform emissivity e t since evapo- 
rated metal surfaces have high (90-99 percent) reflectivity (Reference 3) and since only 
a few reflections are involved. Figure 2 shows the field of view for various thermistor- 
to-cone radius ratios r l /r 2 and aperture angled. 

The energy exchange is judged better by a second mode of imaging, that of imaging 
the source (Figure 3). The detector surface appears completely surrounded by the tar- 
get. Radiation exchange between thermistor and earth is then described by 

P = e t ^(T e 4 - T t 4 ) a t . (2) 

The term = 2ttt^ (I - cos 0 ) represents the detector area. Without the conical 
reflector the same detector would receive approximately P sin 2 ^. The gain introduced 
by the cone becomes about 5 for 0 = 25 degrees. The same method of imaging permits 
estimation of the amount of radiation between cone and detector, as shown in Figure 4. 
Since the emissivity of the cone material can be kept as low as 0.02 and since a large 
portion of the cone will be at the satellite temperature, this effect — small to begin 
with — can be incorporated into the radiation exchange of the unexposed side of the 
thermistor mount. 
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Figure 4 - Radiation exchange between cone and detector 


RESOLUTION AND BANDWIDTH 


The resolution of such a system and the electrical bandwidth required to transmit 
the gathered information are related to the altitude and speed of the vehicle. The de- 
tector, slowly scanning the earth as the satellite moves along, responds to the average 
radiation within its field of view (Figure 5). Any arbitrary intensity distribution in the 
scanning direction can be taken into consideration by superposition of periodic distribu- 
tions of wavelength \ n : 


I, = I*,, eX P • (3) 

The center x 0 of the field of view travels with subsatellite speed v 0 as the satellite, 
visualized as a stable platform, orbits around the earth. Consequently, x 0 is equal to 
vt; and, since 

X = X Q + COS 0 f 

the radiation received by the detector of radius from a solid angle n is given by 

P n = J r 1 2 a n exp (ik n vt) exp (ik n A, cos 0) dH , (4) 

where k n stands for 2n/k n . The increment of the solid angle dfi is derived from geo- 
metrical considerations: 
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dfi 


cos 


x d<fi 

h 2 + ' 


(5) 


Equation 4 can be integrated over 4 > by Sommerfeld's integral for Bessel functions (Ref- 
erence 4): 


>.-r 


2ttt * a n exp(ikvt) h J 0 (k / i)(h 2 + X^) x l dx 


(6) 


This integral can be solved by partial integration: 


P n “ 111 \ a n ex P(i^ v ^) ZZ^TT 

m = 1 q- 1 


(2q - 1) cos 0 sin 2m 8 • 


(7) 


The first null in Equation 7 characterizes the resolving limit of the system. A small, 
square field of view exhibits a sin x/X characteristic and shows the first null at X x = tt. 
A small, circular field of view has a 2J 1 (X)/X scanning function, and the first null is at 
X x = 3.8. Evaluation of the sum in Equation 7 yields X : = 4; an altitude of 300 miles and 
an angle 0 of 25 degrees were assumed. The wavelength X , which corresponds to the 
first null, is then 


2? jH 


( 8 ) 


The corresponding electrical frequency band is given by Equation 9 (when v = 4.4 miles 
per second, and R = 140 miles): 



Figure 5 - Detector slowly scans earth as satellite moves along vector v, 
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The thermal time constant of the detector has to be chosen to accommodate to the highest 
frequency. Compared with those of other detection systems, these bandwidth require- 
ments are very low. 


ENERGY TRANSFER WITHIN THE CONE 


A cross section of the detector assembly (Figure 6) shows the thermistor imbedded 
in the central area on the tip of the cone. The thermistor material could be deposited 
directly on the lower inner part of the Mylar cone. For small angles 0 the black de- 
tector, being truly gray since its emissivity is about 0.8, would resemble an ideal black- 
body. This procedure is not too advisable for the white detector, since the reflectivity 
for visible light is also reduced. The cover plate that shields the cone or cones from the 
interior of the satellite is made highly reflective (e = 0.02). However, in the calculations 
an emissivity of 1 (e s =1) was assumed. This seems to be justified for the following 
reasons. 

Because of multi-reflections, the effective emissivity of two parallel surfaces of 
emissivity e x and e 2 facing each other is 


If both emissivities e x and e 2 are small and equal, e eff becomes e/2. However, if 
radiation emitted by the rear surface of the thermistor is not reflected back to the de- 
tector but bounces back and forth between cone and cover plate until finally absorbed, 
then the cavity between cone and cover acts like a blackbody independent of wall emis- 
sivity, and e eff becomes equal to e rather than to e /2- The actual physical structure 
will be somewhere between these two cases, probably closer to the black cave model. 
Therefore, for the purpose of this calculation, and also to incorporate a small amount of 
radiation from the front surface of the detector to the outer areas of the cone, an emis- 
sivity of 1 was chosen for the cover plate. 

The energy flux Q(r) in the cone material is proportional to the temperature 
gradient, 


Q(r) 


- cD 277 r sin 6 


dT(y) 

dr 


(ii) 


Simultaneously, Q(r) is the flux into the cone tip plus the received or emitted radiation: 



EARTH 



Figure 6 - Energy transfer within conical detector 
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Q(r) = Q( r , ) + 


j eCT ( T e 4 ' T y) da " \ ^( T y 4 ~ T s 4 ) da 


(12) 


The balance between incoming and outgoing radiation of the thermistor is 


Q(rj) = € t a( T e 4 - T t 4 ) 2rrr ? (1 - cos 6) 


, , - T t) 

ea(T 4 -T 4 )da J • (13) 


These equations can be solved, after linearization by means of a reference temperature 


Ty = T o + y 


T y 4 = T 4 + 4T 0 3 : 


T t = T 0 + t 


T 4 = T 4 + 4T„ 3 t + 


T s = T 0 + S 


T s 4 = T 0 4 + 4T 0 3 S + 


T e = T 0 + E 


T e 4 = T 0 4 + 4T 0 3 E + 


(14) 


The error introduced by the omission of second-order terms can be kept small if T 0 is 
chosen close to the mid-range value of the temperatures involved. Furthermore, with- 
out loss in generality, the satellite temperature s can be taken equal to unity and the 
equivalent earth temperature equal to minus unity. Equations 11 through 13 then become 


Q(r) = - cD 2 t 7 r sin 6 


dy 

dr 


(15) 


Q(r) = Q( r x ) - 2eo27T sin 8 4T 0 I y r dr 


• f . 


(16) 


Q( r x ) = 


- e t 2 (1 - cos 0)(1 + t) + e (1 - t) sii 


B + 


2c_,q„(l - t) 


• (17) 


Differentiation of Equations 15 and 16 leads to a Bessel differential equation: 


y + — y 


2e<r4T 0 3 

cD 


y - o 


(18) 


which is solved by Bessel and Hankel functions of order zero and imaginary argument 
(Reference 5), 


y = aj 0 (i£r) + bH^ X) (ifr). 


(19) 
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The parameter £ is introduced for brevity: 




2 


2 


<e g T 0 3 
cD 


( 20 ) 


Boundary conditions determine the constants a and b in Equation 19. One boundary 
condition states that the temperature of the base of the cone must be the satellite 
temperature: 


y(r = r 2 ) = +1 . (21) 

The introduction of Equation 21 into Equation 19 yields 

_ H 0 (i£r) a[j 0 (ifr)H 02 - J 02 H 0 (igr)] 
y ' «o 2 H 02 • (22) 

At the other boundary r = rj, the thermal flux must be Q(r 2 ) as specified by Equation 
17, which determines the remaining constant a : 


(fr 


i 


4e sin 6 


e(yO + sin 0) ll - 


^01 
^0 2 


- e t 2(l - cos 0) 


' “0 2 ' j “02 


4e s 


^Yn~0 [ e(p + sin + £ t2d - cos 5)j (J q 1 H 02 - J 02 H 01 ) + i(J n H 0J - J oa H 01 l 


. (23) 


In Equation 23, p stands for 


p 


C w^w 

ieor 4 T 0 3 ti r j 2 


(24) 


Since the general solution for y equals t for r = r x , the thermistor temperature is 
given by 


»01 + a ^Jpl H 0 2 ~ J02 H 0l) 

H 02 ^02 


(25) 


and finally, by substitution of a and rearrangement of terms, 


t 


e(p + sin 6) - e t 2(l - cos 6) J n H 01 - J 01 H n 

4e sin 6 ^ r i Jqi h 02 I Jq 2 H o i 

e(p + sin 0) + e t 2fl - cos 0 ) A JnH Q2 - J 02 H n 

4e Sin ^ ^ r l Joi H 02 “ Jq2 H 01 


(26) 


Figure 7 shows the thermistor temperature t expressed by this equation (using numeri- 
cal values from Reference 4). The thermistor temperature is within the two limits: the 
satellite temperature S = + i, and the earth temperature E = - 1. A variety of parame- 
ters consolidated in x x and x 2 determines t . The abscissa in Figure 7 is x x , and the 
curves correspond to a particular value of x 2 . The parameters x x and x 2 are 



= 100 



and 
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proportional to the dimensions of the thermistor, r 19 and of the cone, r 2 . They should 
be chosen so that t approaches the earth temperature E as closely as possible without 
unnecessarily increasing the physical dimensions of the detector. The value Xj = 3 
seems to be a good compromise. For thermal reasons alone, x 2 > 6 would be sufficient; 
however, a well defined field of view requires * 2 > 4x r The radii t x and r 2 can then 
be determined: 

r ‘ = Xi ! = * ! vS? ' (27) 

Numerical values for a 1-mil gold-coated Mylar cone yield a thermistor size of about 1 
square centimeter, a reasonable dimension. 

To this point, the influence of wires connecting the thermistor electrically to the 
resistance-measuring instrument has been neglected. This is justified, as can be seen 
from Equation 26, as long as p < sin 0 or c w q w < ^4T 0 3 7rr 1 2 sin 6 , a condition that can 
easily be met by a thin platinum wire (gage number 30, l - r 2 ).. Heat conduction through 
the wires can therefore be neglected. 


CALIBRATION OF THE INSTRUMENT 


A very important task will be the determination of the coefficients k . The linearized 
form of Equation 1 for both the black and the white thermistors yields two equations: 

t bl = k ll E + k 12 S + k 13 A ’ 

t w = k 21 E + k 22 S + k 23 A. (28) 

The constants k n , k 13 , and k 21 will be chosen high, perhaps 0.9. By proper coating, 
the residual responsivity of the white detector to sunlight k 23 will be made as small as 
possible. Calculations in this paper have shown that the combined effect of radiation and 
conduction from the satellite structure, k 12 and k 22 , can be made very small indeed 
(0.1). The experimental verification and determination of all k's is the most important 
part of the calibration. The temperatures of the black and the white thermistors and the 
satellite structure will be telemetered back to earth; then Equation 28 can be solved for 
the unknowns: the blackbody temperature E and the reflected solar energy A. 
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Appendix A 


LIST OF SYMBOLS 


A 

a t 

a, b, a, k 
c 
c w 

D 

da 

E 

f 


( 1 ) 

(ifr) 


*01 


l 02 


H 


11 


h 

I 


albedo 

detector area 
constants 

specific conductivity 
specific conductivity of connecting wire 
thickness of cone material 
areal element 

temperature of earth relative to reference temperature 
frequency 


Hankers function of first kind, m th order of the argument ( igr ) 


Ko'Wl) 

H rWa> 


H j 1 } (if r l) 


satellite altitude 
intensity 

Bessel’s function of m th order of the argument (ifr) 


Joi Jof^ r i> 

Jo 2 Jo^ i( ^ r 2^ 

Jll Ji(l£ r i) 


277 

c - 

•e 

P 

Q(r) 

% 
q, m, n 


wave number 

length of connecting wire 
total radiant flux 
thermal flux in cone 
cross section of connecting wire 
indices of summation 
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R 


r i 

r 2 
r 

S 

k 

T 

T e 

T 0 

T s 

T t 

T y 

t 

t 

v 

X 1 

X 


X 1 
x 2 
y 

P 


n 

p 

a 

e 

€ 


radius of target area 

radius to arbitrary point in the target area 
radius of thermistor 
radius of cone 

radius to an arbitrary point of the cone 

temperature of structure of satellite relative to reference temperature 

solar constant 

absolute temperature (°K) 

equivalent temperature of earth 

reference temperature 

temperature of mounting structure 

temperature of thermistor 

temperature of cone material 

temperature of thermistor relative to reference temperature 
time 

subsatellite speed 
first zero F(X x ) = 0 
subsatellite path; arbitrary variable 
center of field of view 
auxiliary parameter (see Figure 7) 
auxiliary parameter (see Figure 7) 
temperature of cone relative to reference temperature 
elevation angle 

emissivity of cone surface (= 0.02) 
emissivity of mounting structure in satellite (= 1) 
emissivity of thermistor (= 0.8) 
wavelength of intensity distribution (Figure 5) 
solid angle 
azimuth angle 

c qyi ecr4T 0 3 nr * as defined by Equation 24 
Stefan- Boltzmann radiation constant 
aperture angle (one-half of field of view) 
y2ea4T 0 3 /cD as defined by Equation 20 
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